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ABSTRACT 

The thermal Sunyaev-Zel'dovich (tSZ) effect directly measures the thermal pressure of free electrons inte- 
grated along the line of sight and thus contains valuable information on the thermal history of the universe. 
However, the redshift information is entangled in the projection along the line of sight. This projection effect 
severely degrades the power of the tSZ effect to reconstruct the thermal history. However, this otherwise lost 
redshift information can be recovered by the SZ tomography technique, namely by cross correlating the tSZ ef- 
fect with galaxies of known redshifts, or alternatively with dark matter distribution reconstructed from lensing 
tomography. For this purpose, we investigate in detail the 3D distribution of the gas thermal pressure and its re- 
lation with the matter distribution. We compare our adiabatic hydrodynamic simulation with the one including 
radiative cooling and star formation with supernova (SN) feedback. We confirm that these additional processes 
of gastrophysics suppress the pressure power spectrum and thus the SZ angular power spectrum significantly. 
On the other hand, they only affect the 3D pressure-matter correlation coefficient r(k) at ~ 1 — 2% level. Since 
r(k) is a key input in the SZ tomography, this finding makes the tSZ tomography robust against uncertainties in 
the understanding of cooling and SN feedback. 

Subject headings: cosmic microwave background - large-scale structure of universe - galaxies: clusters: general 
- methods: numerical 



1. INTRODUCTION 

Free electrons in the universe reveal their existence in 
the CMB sky through their inverse Compton scattering of 
CMB photons. The induced secondary CMB temperature 
anisotropies, proportional to the electron thermal energy in- 
tegrated along the line of sigh t, are the well known ther- 
mal S unyaev- Zel'dovich effect (ISunyae v & Zeld ovichlll972L 
Since massive galaxy clusters contain large reservoir 
of hot electrons, the generated tSZ effect can thus overwhelm 
the primary CMB around cluster scales. For this reason, 
the thermal SZ effect of dozens of galaxy c lusters has been 



meas u red by various experiments (refer to iCarlstrom et al 
2002 1 : iReese et al.1 120021: Uones et al.l 120051: lLaRoque etal 
2006t IBonamente et al.ll2006l for reviews). On the other hand 



blindly detecting the tSZ effect in random directions of sky is 
much more difficult, since the expected signal is overwhelmed 
by the primary CMB fluctuations. Currently there are only 
a few tentative dete ctions through the observed small scale 
CMB power excess (iDawson et ai]l2002L 120061: iMason et al 
Goldstein et al] 120031: iRunvan et al.l 120031: iKuo et al 
2004t|Readhead et al . 2004; B ond et al.ll2005HReichardt et al 



20080. However, the situation will be significantly improved 



in the next few years by ongoing and proposed ground surveys 
such as SZAEL ACfl APEXQ, SPT0 and the Planck satellite. 

Precision mapping of the SZ sky is of great importance to 
both cosmology and astrophysics. The SZ effect is a power- 
ful finder of galaxy clusters at high redshifts. The efficiency 
of free electrons to generate the thermal SZ effect is redshift 



1 However, refer to the recent SZA measurement for the opposite conclu- 
sion ISharp et al.ll2009l) 

2 http://astro.uchicago.edu/sza/ 

3 http://www.physics.princeton.edu/act/ 

4 http://bolo.berkeley.edu/apexsz/ 

5 http://pole.uchicago.edu/ 



independent. Photons originated from redshift z suffer a fac- 
tor of 1 +z energy loss caused by the cosmic expansion. On 
the other hand, CMB photons that electrons scatter at redshift 
z are a factor of 1 + z more energetic than CMB photons to- 
day. These two effects cancel out exactly and enables galaxy 
clusters to be detected at high redshift without extra effort, op- 
posite to the X-ray cluster finding. An exciting advance in this 
area is the recent discovery of 3 new galaxy clusters thro ugh 
the tSZ effect by the SPT group (IStaniszewski et alJl2008l) . 

The tSZ effect is also a powerful probe to the thermal his- 
tory of the universe, since it directly probes the thermal en- 
ergy of intergalactic medium (IGM) and intracluster medium 
(ICM). Numerical simulations show that, the amplitude of 
the SZ signal is sensitiv e to the amount of r adiative cool- 
ing and energy fee dback (ISpringel et alj|200it Ida Silva et alJ 
1200 U IWhite et al.|[2002t iLin et aUl2004l) . However, it is not 
straightforward to extract information on these astrophysical 
processes from the tSZ measurements alone. First of all, re- 
search shows that there exist great degeneracies between dif- 
ferent competing processes. Even worse, the tSZ effect only 
measures the electron thermal energy projected along the line 
of sight. The redshift information of these astrophysical pro- 
c esses is thus entangle d in the projection. 

IZhang & Pen! d2001l) proposed to recover the redshift infor- 
mation by cross correlating the tSZ effect with galaxies with 
at least photo-z information. The idea is that galaxies in a 
given redshift bin should strongly correlate with the tSZ signal 
from the same redshift bin. A key link between the measured 
cross correlation and the gas pressure auto-correlation that we 
want to extract is the cross correlation coefficient r between 
the thermal energy and the galaxy number density. Assum- 
ing a constant r, the time resolved thermal energy distribution 
can be reconstructed self consistently. This SZ tomography 
technique would be applicable in reality, since SZ surveys of- 
ten have follow-ups of galaxy surveys. For example, the dark 
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energy surve}0 will cover the SPT sky and measure photo- 
metric redshifts of <~ 10 8 galaxies up to z = 1.3. The lensing 
tomography also helps to reconstruct the 3D matter distribu- 
tion, which can also be correlated with the SZ map to make 
the SZ tomography. 

In the current paper, we reformulate this SZ tomography 
technique and explore the possibility to improve its robust- 
ness. We no longer approximate r as a constant. Rather, 
we rely on numerical simulations to quantify its scale and 
redshift dependence. We are able to show that, r is insensi- 
tive to gastrophysics such as radiative cooling and supernova 
(SN) feedback. Namely, r calculated from adiabatic hydrody- 
namic simulations should be sufficiently accurate, even with 
the presence of radiative cooling and SN feedback. We are 
then able to take this r as input to perform the SZ tomogra- 
phy. 

The paper is organized as follows. In §J2] we introduce 
the tSZ effect and explain the SZ tomography technique. We 
then analyze our high precision hydrodynamic simulations to 
quantify the dependence of various SZ statistics on the extra 
gastrophysics in $3] Although most quantities are sensitive to 
these gastrophysical processes, we find that the cr oss c orrela- 
tion coefficient r only weakly depends on them ( §13.31 ). This 
feature allows us to use r calculated from the adiabatic simula- 
tions as the input of the SZ tomography, despite the existence 
of complicated gastrophysics in the real universe. We discuss 
and make conclusion in ^4] 

2. THE THERMAL SZ EFFECT AND THE SZ TOMOGRAPHY 

The tSZ effect induces a new source of CMB temperature 
fluctuations with the amplitude 

Ar(0) 



TcMB 



■g(x)y(0). 



(l) 



Here, 9 is the direction on the sky. g(x) describes the spectral 
dependence. In the non-relativistic limit, 



g(x) = 



e*+\ 
V-l 



(2) 



where x = hv/k B TcMB = ^/56.84 GHz and v is the observed 
frequency of CMB photons. The Comptonization parameter 

m is 

(Tj f 

y= ~ adxn e k B T e , (3) 

m e c L J 

where n e kBT e is the hot electron pressure. \, n e , k B , T e and oj 
are the comoving diameter distance, number density of free 
electrons, the Boltzmann constant, electron temperature and 
the Thompson scattering cross section respectively. 

The spectral dependence of the tSZ effect is unique. The 
tSZ effect shows as CMB temperature decrements at v < 218 
Ghz and as increments at v > 218 Ghz. The spectral function 
g(x) — > -2 at the Rayleigh- Jeans band i<1 and g(x) — >x— 4 
at x ^> 1. This unique spectral dependence allows a clean 
separation of the tSZ effect from other CMB components in 
multi-band CMB surveys. 

The y parameter contains key information on the thermal 
history of the universe. However, since it only measures the 
projected electron thermal energy along the line of sight, the 
redshift information is smeared by this projection effect. Our 
SZ tomography technique aims to recover the otherwise lost 
redshift information in the tSZ effect. 

6 http://www.darkenergysurvey.org/ 



2.1. The SZ tomography 

One of the most widely used statistical quantities of the 
tSZ effect is the angular power spectrum Cf z . Throughout 
this paper, unless specified, we will focus on the Rayleigh- 
Jeans limit AT /T = —2y. Under the Limber's approxima- 
tion (iLimberll 19541) . Cf z is related to the 3D thermal pressure 
power spectrum Ap(k,z) by the following relation, 



;2 rXcMB 1 

fc; sz =/ AKk=L 



z)W t l z (z) X d X 



(4) 



We have adopted the flat cosmology in the above expression. 
The weighting function is 



Wtsz(z) = -2a T a 



[n e k B T e ] 



m e c e - 



(5) 



The pressure power spectrum A 2 p is that of the fractional ther- 
mal pressure fluctuations 5 P = n e k B T e / (n e k B T e ) - 1. W t $z(z) 
tells us the overall thermal energy of the universe and the A 2 p 
tells us the clustering of the thermal energy. 

Our SZ tomography technique aims to reconstruct the time 
resolved Ap(k,z)W^ z (z). This quantity tells us the overall am- 
plitude of the thermal energy and the clustering strength at 
redshift z. The tSZ effect is correlated with tracers of the large 
scale structure. Given the redshift information of these trac- 
ers, such as galaxies, we can manage to recover the redshift 
information of the tSZ effect. The original tomography is pre- 
sented in the variation formalism dZhang & Penll2001[) . while 
in the current paper, we reformulate it in a more straightfor- 
ward manner. The key idea of the SZ tomography is that, 
galaxies distributed in a certain redshift range correlate with 
the tSZ signals contributed by the IGM in the same redshift 
range. Given the photometric redshift information, we are 
able to split galaxies into different redshift bins. The galaxy 
number distribution in the i-th redshift bin is rii(z), which is 
related to the photo-z distribution nf(z p ) by 



J+Az, 

ni(z)= I p(z\z P )n p (z P )dz f 

-A;,-/2 



(6) 



Here, p(z\z p ) is the photo-z probability distribution function. 
Due to non-negligible photo-z error, the real galaxy distribu- 
tion is wider than Az;. However, it can still be narrow, cen- 
tered at Zi- We normalize n, such that j Q rti(z)dz = 1. The 
tSZ-galaxy cross power spectrum is then 



;2 

' ^tSZ-g 

2^' 



'/ 



A 2 P ^k=-,z)W t sz(z)n i (z) X dz 

X 



(7) 



Notice that the integral in Eq. [7] is over the entire redshift 
range. For a sufficiently narrow redshift bin, the galaxy distri- 
bution n,(z) peaks at z = zQ with a narrow r.m.s. redshift width 
~ y/ ' a 1 + (Azi) 2 1 12, where a z is the photo-z error. For typical 
values a z = 0.05(1 +z) and Az,- = 0.2, the effective r.m.s. width 
is <~ 0.1-0.2. Over this redshift width, the functions Ap ? , W 



7 This relation assumes no catastrophic error in photo-z measurement. 
With the existence of catastrophic error, the peak redshift can be shifted by 
outstanding outliers. 
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and x var Y slowly. We thus have an approximation 



f z - s ^A 2 P (k=-,z i )Wtsz{zdXi 
g Xi 



A 2 Pg (k=-,Zi)W tSZ (zdXi 

Xi 



rii(z)dz 



(8) 



Here, Xi = xfc)- It is clear from the above equation that the 
cross correlation between the thermal SZ maps with galaxy 
distribution picks up and, to a good approximation, only picks 
up relevant information within the given redshift bin. This is 
a key step to recover the redshift information of the thermal 
SZ effect. 

For similar argument, the galaxy angular power spectrum 

is 



I 2 

—C s 
2n l 



dz 



A z (k=-,z)nf(z)x^dz 



o X 
A 2 (k=-, Zi )Xi 



Xi 



' 2/ dz J 

n i( z)-dz. 



(9) 



The next step is to recover the thermal SZ contribution in 
the given redshift bin with the correct weighting, namely, 
Ap(k,z)W^ z (z). This step requires a key input, namely the 
cross correlation coefficient r between fluctuations in the gas 
pressure and the galaxy distribution, defined by 



r(k,z)-- 



Al g (k,z) 



A P (k,z)A g (k,z) 
A 2 g (k,z)W tSZ (z) 



(10) 



[A P (k,z)W t sz(z)]A g (k,z) 
Combining Eq. [8] and Eq. [9] we obtain 
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In this relation, the angular power spectra C, and Cf 
are evaluated at I = kxi- Both the angular power spectra can 
be measured directly by combining the SZ surveys and the 
galaxy surveys. The quantity «,(z) is given by the galaxy sur- 
veys, and x can be calculated given the cosmology model. As 
long as r(k,z) can be figured out, one can apply this equation 
to do the thermal SZ tomography. 

We summarize the thermal SZ tomography as follows. 

• Measure the cross power spectrum cf Z E between the 
thermal SZ effect and the galaxy distribution in a given 
redshift bin. 

• Measure the galaxy power spectrum Cf in this redshift 
bin. 

• Reconstruct the thermal SZ contribution from this red- 
shift bin through Eq. QT| 

A key input requ i red in this SZ tomography is the quan- 
tity r. [Zhang & Pen ( 2001 ) adopted a simplification that r is a 
constant. They further showed that, the value of r can be mea- 
sured combining EqQTjand the measured SZ power spectrum. 



This approach is self consistent and able to provide a quick 
realization of the SZ tomography. However, improvements 
must be made to do precision SZ tomography. As shown in 
33.31 the assumption that r is a constant is only accurate at 
the level of ~ 20%. Furthermore, to obtain the value of r un- 
biasly from the SZ power spectrum, the galaxy surveys must 
well cover the SZ redshift range (z < 2). This is challenging 
for galaxy surveys. 

Since r(k,z) is such a key quantity in the SZ tomography, 
a natural question arises: Can one robustly predict r? This 
seems challenging, given the fact that various complicated 
gastrophysical processes affect the tSZ effect. Surprisingly, 
robust prediction on r is likely feasible. The key point is that, 
these gastrophysics affects Ap and A 2 Pl( in basically the same 
way, so that these effects roughly cancel out in r by defini- 
tion. To quantify the dependence of r on these gastrophysics, 
we compare two sets of hydrodynamic simulations with and 
without radiative cooling, star formation and SN feedback. 
These simulations confirm the above naive speculation and 
find that, these additional gastrophysics can only affect r at 
1% level, despite the fact that they alter the SZ power spec- 
trum by ~ 40%. 

3. SIMULATIONS 

The tSZ statistics has b e en st u died by both semi- 
analytical models JCooravl 120001 IZhang & Pen! 120011: 



Komatsu & Sellakl l200l IZhang & Wul 12003b IZhang et al.l 
.20041 IZhang & Shethl 120071) and numerical simulations 
dPersi et al.ll 19951 IRefregier et al.ll2000 t IRefregier & Tevssier 



20021 Ida Silva et all 120001 120011 iSeliak et al 



Springe! et all 1200 U IWhite et al.l 120021 Lin et all 1200 



Zhang jet all 120041 IRoncarelli et al.l 120071 lHallman et al.1 
20071 120091) . Beyond the gravitatio n al heating mechanism , 
some works (Ida Silva et al l 120011 ISpringel et all 12001 

120041 iRon 



IRoncarelli et al.l 12007 



White et al.l 12002 ; iLin et al 

Scannapieco et al. 120081) have incorporated additional gas- 
trophysics such as radiative cooling, preheating, SN/AGN 
feedback. These studies found that these processes suppress 
the small scale SZ power spectrum. But none of them 
addressed the dependence of the cross correlation coefficient 
r on these processes, whose investigation is the key goal of 
this paper. 

In this paper, we analyze a controlled set of hydrodynamic 
simulations for the relevant SZ statistics. The cosmology 
adopted is a ACDM cosmology: tl A = 0.732, tt = 0.268, Q, b = 

0.04448, A = 0.71 o-g = 0.85 . The simulations are run by 

GADGET2 code (ISpringel I l200l . with 100 h~ x Mpc box 
size, and 512 3 particles for both dark matter and SPH par- 
ticles. The masses of the dark matter and SPH particles are 
4.62 x 10 8 /?"' M Q and 9.20 x 10 7 /?"' M Q respectively. One 
simulation only includes gravitational heating and the other 
includes radiative cooling and star formation with SN feed- 
back. We refer these two as adiabatic and non-adiabatic run, 
respectively. For the non-adiabatic run, 7.63% of gas parti- 
cles turn into stars by z = 0, resulting in stellar mass density 
O* = 0.0034. Both simulations start from redshift z = 120. We 
o utput 60 snapsho ts , which are equal ly spaced in lna. Refer 
to lLinetall (120061) : Uing et al.1 ( 120061) for more details of the 
simulations. 

We derive the temperature 7) of the r-th gas particle as fol- 
lows 

k B Ti = (j-l)nm p iii, (12) 
with fx = 0.588 the mean molecule weight, 7 = 5/3 the spe- 
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FIG. 1 . — The redshift evolution of the density-weighted temperature T p in 
the adiabatic (solid line) and the non-adiabatic simulation (dashed line) with 
cooling, star formation and SN feedback. 



cific heat ratio for monatomic gas, and m, the internal energy 
per unit mass of the 2-th particle. When the temperature of a 
SPH particle is higher than 10000 K, we consider it as ion- 
ized, otherwise as neutral. We sum up all the SPH particles to 
calculate the density-weighted temperature T p at all the avail- 
able snapshots. From these mea surements, we are then able 
to calculate the y parameter §13.11 

We use the fast Fourier transform (FFT) to calculate the 
power spectra of dark matter, gas pressure and the cross power 
spectrum between the two. For dark matter particles, we as- 
sign the masses to 512 3 grids using the Clouds-In-Cells (CIC) 
algorithm. On the other hand, we assign the thermal energy of 
SPH particles to the same grids using the SPH kernel while as- 
suring the energy conservation. The gas pressure on the grid 
is converted from the gas thermal energy assigned. We re- 
move the shot noises from the corresponding raw power spec- 
tra. However, since we are mainly interested in scales well 
below the Nyquist frequency (e.g. k < 16 h Mpc -1 ), we don't 
correct for the alias effect^ Results are shown in {s!3.2l & [3~3l 

3.1. The tSZ mean y parameter and the gas density weighted 
temperature 

The mean Comptonization y parameter is given by the fol- 
lowing equation 

a j f _ 

y= r / n e k B T p adx- (13) 

m e c- J 

8 The alias effect would affect the power spectra to scales as large as a 
quarter of t he Nyquist frequency, nam ely k ~ 4 h Mpc -1 in our simula tions 
<Jing|2005D . The methods proposed bv Oml (2003) and lCui et al] <2008l) can 
effectively correct for this effect. However, to the first order approximation, 
it affects the corresponding power spectra in the two simulations in the same 
manner. Since we're mainly concerned about the differences between the two 
power spectra, we neglect this effect in this paper. 



FIG. 2. — The distribution of internal energy at different temperature u(T) 
which is defined as J u(T)dlnT = U, where U is the total internal energy. The 
additional processes brings the total internal energy down by 12.0%, 10.9%, 
1 1.5% and 16.6% at z=0, 0,52, 1.02 and 2.08 respectively. The same arbitrary 
unit is adopted for both the simulations. For both the adiabatic (solid lines) 
and non-adiabatic (dashed lines) simulation, the most contributions come 
from the gas at temperature around 1-2 keV. And the lower the redshift, 
the more contributions from higher temperature components. 

Here T p = (n e T)/n e is the density weighted temperature, 
which is computed by the following equivalent expression 

Y"" . T:tn: 

T p = ^-^ . (14) 

The sum is over all gas particles, assuming gas with T, > 
10000 K as ionized, otherwise neutral. 

Figure Q] shows the redshift evolution of the density 
weighted temperature for the adiabatic (solid lines) and 
the non-adiabatic (dashed lines) simulation. As shown, 
the additional gastrophysical processes suppress the mass 
w eighted temperature T p , which confirms the earlier results 
bv lWhite et alj d2002l) . though more evident than theirs at z=0. 

We further look into the effect of cooling and feedback in 
redistributing the energy budget into regions with different 
temperatures. The logarithmic contribution to the total ther- 
mal energy from regions with different temperatures is shown 
in Fig. |2] We can see clearly that, for the 4 presented redshifts, 
the most contributions of the total energy come from the re- 
gions with T ~ 1 - 2 keV at lower redshifts except z = 2.08, 
for both adiabatic and non-adiabatic simulations. The figure 
also shows that gas with temperature higher than ~ 2 keV 
is slightly affected by these additional gastrophysics. Most of 
these gas lies in clusters, which have deep gravitational poten- 
tial wells to offset the feedback effect. Also, the star formation 
rate is relatively low in cluster regions, resulting in less feed- 
back intensity and also less depletion of hot gas by cooling. 
On the other hand, gas cooler than ~ 1 keV has opposite fate. 
Most of these gas lies in galaxy groups or in IGM , susceptible 
to cooling and feedback. This suppresses their contribution to 
the total energy. The combined effect is that, the total thermal 
energy is suppressed by a factor of 12.0%, 10.9%, 1 1.5% and 
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FIG. 3. — Top panel: The power spectra of gas pressure Ap(k) at red- 
shifts z=0,0.52,1.02 and 2.08. At small scales, the pressure power spectra 
of non-adiabatic simulation (dashed lines) are greatly reduced by cooling 
and SN feedback, even more significant to higher redshifts, compared to the 
adiabatic one (solid lines). However, the powers at large scales are slightly 
boosted. Bottom panel: The differences between the power spectra defined 
by A(A 2 p )/Al A = (Al NA (k)-Al A (k))/A 2 pA (k) 

16.6% at z=0, 0,52, 1.02 and 2.08 respectively. 

Given T p measured at the 60 simulation snapshots, we are 
able to integrate it over the whole relevant redshift range 
to calculate y, according to Eq. [13] We adopt a simple 
linear, while stable, interpolation to model T p at redshift 
falling between two adjacent output redshifts. For both cases, 
about a half of the y signal comes from z < 1. We obtain 
y = 2.07 x 10" 6 for the adiabatic run and y = 1.77 x 10" 6 
for the non-adiabati c run. These values of y are lower than 
the results in e.g . (Ida Silva et al.l 1200 U IWhite et al.l 120021: 
IZhang et al.ll2004l) . so as the result of T p . On the other hand, 
it is higher than some rec ent simulation results, such as that 
of lRoncarelli etal.1 ( 120071) using GADGET2. Differences in 
cosmological parameters, such as difference in a% and Qj,, 
can only account for part of the discrepancy, while cosmic 
variance, differences in the input gastrophysics and numeri- 
cal codes may be responsible for the remains. We postpone 
further investigation of this issue elsewhere^ 

3.2. The 3D gas pressure power spectrum and the tSZ power 

spectrum 

Most of the simulated SZ power spectra in the literature 
are directly calculated from the simulated 2D SZ maps. In 
this paper we take a n alternative approach, as a dopted by 
IRefregier et"aT1(l2000l) : lRefregier & Teyssierl (120021) . For each 

9 We thank Volker Springel for providing the result of an independent 
GADGET2 simulation, which allows us to do a preliminary check. The cos- 
mological parameters of this simulation are identical to ours, except the dif- 
ferent initial conditions and ag = 0.9. This simulation gives T p = 0.283 keV. 

Taking the scaling relation T p tx o"g 05 jZhang et al.ll2004l) . we esti- 
mate that this simulation would result in T p = 0.243 keV for ag = 0.85, in 
good agreement with ours. 
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FIG. 4. — The bias bp(k) = y Ap(k) / A} n (k) of gas pressure relative to the 
dark matter at the four redshifts. The same legends are adopted as in Figure 

m 

simulation output, we can calculate the 3D power spectrum of 
the gas thermal pressure Aj>, The 2D SZ power spectrum is 
obtained through Eq. d4]i and we interpolate linearly between 
adjacent snapshots to model A^(fe) at any other redshifts. 

A£ at z = 0, 0.52, 1.02 and 2.08 are shown in Figure[3] At 
small scales, the pressure power spectra of the non-adiabatic 
simulation (dashed lines) fall below the corresponding adi- 
abatic ones (solid lines). Most of tSZ signals at these scales 
are contributed by gas in less massive halos. These halos have 
shallower potential wells and thus SN feedback is easier to de- 
plete gas. This suppression effect is even larger toward higher 
redshifts where the gravitational potential wells are shallower, 
leading to as high as 50% reduction at z = 2.08. The same 
feedback redistributes gas and moves the clustering power of 
the thermal energy from small scales to large scales. This 
is likely the reason that we see enhancement at large scales 
(Fig. [5]). However, at these scales, massive halos contribute 
most of the signal. Since they have deeper potential wells to 
fight against feedback and confine gas, the large scale pressure 
power spectrum is less affected. 

The gas pressure power spectrum is often expressed by the 
gas pressure bias with respect to the matter distribution, de- 
fined as bp(k) = AUk)/ Al t (k). The measured pressure bias 
is shown in Figure [4] We find in Figure H] that the results 
are in good agreement with the simulation results at z = by 
IRefregier & Teyssierl (120021) . 

The SZ power spectrum evaluated from the Limber integral 
is shown in Figure In the existence of the cooling and self- 
regulated star formation, the tSZ angular power spectrum at 
large scales is suppressed by ~ 20%, which is mainly due to 
the reduction of density weighted temperature T„, At the same 
time, the powers are reduced by as large as around ~ 40% at 
very small scales around / = 10000, which is jointly caused 
by the suppression of both the density weighted temperature 
and pressure power spectra compared to the adiabatic sim- 
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FIG. 5. — The SZ angular power spectra C ( tSZ for the two simulations. Ra- 
diative cooling and SN feedback suppress the SZ power spectrum by ~ 20% 
at large scales and around ~ 40% at small scales. The same legends are 
adopted as in Figure [T] 

ulation. These suppression effects confirm the previous re- 
sults by|dI~Silvi~iLS(|20^ And what's 
more, the non-adiabatic power spectrum peaks on a slightly 
larger scale, about / ~ 8000, than the adiabatic one, which 
may arise from SN feedback that would d rive away the gas out 
to diffuse regions. ISpringel et al.l (1200 lh also found that en- 
ergy injection would suppress the po wer on small sca l es and 
push the powers toward larger scales. iRoncarelli et al.l (12007b 
considered the same non-gravitational processes as our non- 
adiabatic one in their simulation, and their power spectrum 
behaves consistently with ours, despite of their low y. 

3.3. The correlation coefficients r(k,z) 

As addressed in $2] a key input of the SZ tomography is 
the gas pressure-galaxy number density cross correlation co- 
efficient r. State of art numerical simulations with purely 
gravitational heating (adiabatic simulations) is able to model 
r(k,z) robustly. If we can further quantify the dependence of 
r(k,z) on additional gastrophysics such as radiative cooling 
and feedback, we will be able to model r(k,z) robustly for 
general cases and perform the SZ tomography robustly. 

For this purpose, we analyze the behavior of r(k, z) in the 
two simulations. Since galaxies are tightly correlated with the 
matter distribution, and the stochasticity is small for the mea- 
surement of r(k, z) between the gas pressure and galaxy num- 
ber density, we are able to approximate galaxy distribution as 
matter distribution. Namely, what we actually measure is rp m , 
the cross correlation coefficient between the gas pressure and 
the matter density. To the first order approximation, these two 
should be identical. For example, for the often adopted con- 
stant galaxy bias model, the two are equal to each other. Even 
for scale dependent galaxy bias, as long as the stochasticity 
is negligible, the two should be equal to each other. In future 
works, we will populate dark matter halos with galaxies and 
measure the real gas pressure-galaxy cross correlation coeffi- 




k/h Mpc 



FIG. 6. — Top panel: The cross correlation coefficient between the 3-D 
gas pressure and dark matter in adiabatic simulations (solid lines ) and non- 
adiabatic simulations (dashed lines). For both the simulations, r(k) is greater 
than 0.9 at k ~ 0.062 h Mpc -1 for all the 4 redshifts, as gas pressure traces 
dark matter on the very large scales. At intermediate scales where 0.4 < k < 3 
h Mpc~' , there's still great correlation about ~ 0.8 at z = 0, while r(k) drops 
by 25% to ~ 0.6 at z = 2.08. Bottom panel: The difference of r between 
the two simulations: Ar(k) / r&(k) = (r N ^(k)-r^{k))/r^(k). If cooling and SN 
feedback are included, there's only ~ 1 - 2% change in r(k) at most scales, 
and only at the smallest scales does r change ~ 5%. 



cient. 

The result on r(k,z) is shown in Figure[6] We find gas pres- 
sure tends to trace the matter distribution faithfully on the very 
large scales. At the largest scale that our simulation can ap- 
proach, k~2ir/L= 0.062 h Mpc" 1 , r(k) > 0.9 for all of the 
4 redshifts. At 0.4 < k < 3 h Mpc" 1 , there's still great cor- 
relation between pressure and matter, while r(k) changes by 
25% from ~ 0.8 at z = to 0.6 at z = 2.08. However, at each 
redshift, the correlation strength nearly keeps unchanged in 
this range. At smaller scales k > 3 h Mpc" 1 , r(k) begins to 
decrease faster. 

We further investigate the influences caused by cooling and 
star formation on r(k). Surprisingly, we find that these addi- 
tional gastrophysical processes have little impact on r(k). In 
most of the relevant k range, the change is ~ 1 %. Only at very 
small scales (corresponding to / > 10 4 ), can the effect reach 
5% level. This behavior is however explainable, since the gas- 
trophysics influences the denominator and numerator of r(k) 
in basically the same way. A 1% error in r results into a 2% 
error in the reconstructed ApW,| z , negligible to the ~ 30% 
change in ApW t | z induced by cooling and feedback. There- 
fore, this result enables us to extract robust predictions for 
r from adiabatic simulations without considering additional 
complicated gastrophysics, and then carry out the tSZ tomog- 
raphy method by applying r back into the observations. Since 
eventually we need a large sample of simulations with high 
resolution to improve the calculation, this feature can signifi- 
cantly save the computation cost. More importantly, it breaks 
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the otherwise circular procedure to do the tSZ tomography^ 
This feature assures us that we can rely on adiabatic simula- 
tions alone to provide r. However, it does not mean that the 
simulations presented in this paper have already robustly pre- 
dicted r. We expect the function r(k) to be sufficiently smooth 
with respect to k. However, since our box size is only lOO/i -1 
Mpc, statistical fluctuations induced by cosmic variance are 
non-negligible and the simulated r-k relation shows clear ir- 
regularities (Fig. [SJ. To beat down the cosmic variance, we 
need more simulations or larger box size. Nonetheless, we 
only need adiabatic simulations for this improvement. 

3.4. The galaxy-matter relation 

The r measured above is actually between the gas pressure 
and the matter distribution. We have shown its robustness 
against various gastrophysical processes. However, in the pro- 
posed SZ tomography, what we really need is the one between 
the gas pressure and the galaxy distribution. At sufficiently 
large scale, the galaxy distribution with respec t to the matter 
distribution is believed to b e deterministic (e.g. iBonoli & Pen! 
120091: lBaldaufeTaT]l2009h . This is a generic behavior in the 
standard model of structure formation and is unlikely to be 
changed by realistic gastrophysics. Thus at these scales, our 
result proves that r between gas pressure and galaxies is also 
robust various gastrophysical processes. However, at small 
scales, stochasticities between galaxy and matter distribution 
develop and the cross correlation coefficient between galaxy 
and matter distribution r mg ^ 1 . In this case, we should also 
quantify the impact of various gastrophysical processes such 
as feedback and cooling on r mg . 

Unfortunately, since even the state of art simulations are 
not able to simulate galaxies from first principles, we are not 
able to robustly quantify the impact of these gastrophysical 
processes on r mg . Instead, we will adopt the halo model and 
test r mg by varying parameters of the halo occupation distribu- 
tion (HOD) to study possible impact of these gastrophysical 
processes. The technical details and results are presented in 
the appendix. We find that in some cases r mg barely changes 
and in other cases it can change by ~ 10%. Due to uncer- 
tainties in the dependence of HOD parameters on these gas- 
trophysical processes, we are not able to draw any conclusive 
observation. However, these results do warn us the possibility 
that these astrophysical processes could significantly change 
r mg and hence the cross correlation coefficient between the 
gas pressure and galaxy distribution, despite the fact that they 
have only minor impact on the cross correlation coefficient 
between the gas pressure and matter distribution. 

The good news is that, even for this worst scenario, we are 
still able to carry out the thermal SZ tomography robustly, 
with the help of the lensing tomography. SZ surveys of- 
ten overlap with lensing surveys. For example, the South 
Pole Telescope survey basically overlaps with the Dark En- 
ergy Survey (DES). Through the lensing tomography of sur- 
veys with redshift information, in principle we are able to 
reconstruct the 3D matter distribution from the nearby uni- 
verse to z ~ 2 whe re the furtherest source galaxies reside. 
iMassev et al.l d2007l) have clearly demonstrated the feasibility 
of this technique. Given the capacity of stage III and IV lens- 
ing surveys, the 3D matter distribution can be reconstructed 

10 If r(k) strongly depends on the detailed gastrophysics, the tSZ tomog- 
raphy procedure will become circular. In this case, we need to know the de- 
tailed gastrophysics to reliably predict r(k, z) from simulations and then apply 
back to observations, while we count on observations to tell us the detailed 
gastrophysics. 



with much higher accuracy. We are then able to cross corre- 
late the thermal SZ measurement with the reconstructed mat- 
ter distribution in a redshift bin to perform the thermal SZ 
tomography. The thermal SZ tomography can then be carried 
out as follows. 

• Reconstruct the 3D matter distribution from the lensing 
tomography. 

• Measure the cross power spectrum Cf z ~ m between the 
thermal SZ effect and the matter distribution in a given 
redshift bin. 

• Measure the matter angular power spectrum Cf in this 
redshift bin. 

• Reconstruct the thermal SZ contribution from this red- 
shift bin through the following relation 

A 2 P Jk, Zi )W t sz 



A' P (k,Zi)W^ z ( Zi ) = r~\k,zd 



-r-\k,Zi) 



A 2 m (k,Zi) 

/ 2 [cf^] 2 / oo nfe)£<fe 

2-KCf X i 



(15) 



In this case the relevant cross correlation coefficient r is the 
one between the gas pressure and the matter density, exactly 
the one that we study in this paper. We have demonstrated its 
robustness against various gastrophysical processes and hence 
the robustness of the thermal SZ tomography. 

4. CONCLUSION AND DISCUSSION 

In this paper, we investigate the feasibility of the tSZ to- 
mography method, which aims to extract the redshift infor- 
mation of the tSZ effect. As explained in $2] future SZ 
and galaxy surveys will measure the galaxy angular power 
spectrum Cf in each redshift bin and the corresponding tSZ- 

galaxies cross power spectrum cf Z g . Given the cross corre- 
lation coefficient r(k,z) between the gas thermal energy and 
the galaxy number density, we are able to reconstruct the 
time resolved thermal energy distribution. So the key ques- 
tion in the tSZ tomography is whether we can robustly pre- 
dict r( k, z), given uncertainties in gastrophysical processes. To 
quantify the effect of two dominant gastrophysical processes, 
namely radiative cooling and SN feedback, we compare the 
result of our adiabatic simulation against the one with cool- 
ing and feedback. We find that r(k,z) is insensitive to these 
additional gastrophysics. The resulting difference in r(k,z) is 
~ 1-2% in most relevant k and z range, much smaller than the 
~ 30% change in the SZ power spectrum and the gas density 
weighted temperature. This allows us to neglect the depen- 
dence of r(k, z) on these gastrophysical processes and adopt 
the r(k,z) evaluated from adiabatic simulations as the input of 
the tSZ tomography. We thus show that the tSZ tomography 
is feasible in reality. 

There are a number of improvements which can be made in 
future work. (1) As addressed in the end of ^3] We need to run 
more simulations to reduce the cosmic variance. We may also 
need to correct for the alias effect to improve the accuracy of 
r at scales k>4h Mpc -1 . (2) In this paper, to evaluate r(k,z), 
we have approximated galaxy distribution as dark matter dis- 
tribution for simplification. In the limit that the stochastic - 
ity vanishes, our approximation becomes exact. This should 
be the case at large scales. At small scales, non-vanishing 
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stochasticity could result in some systematical errors in the 
tSZ tomography. In the future work, we will populate halos 
in our simulations with galaxies and quantify the possible sys- 
tematics induced by this simplification. (3) Another issue is 
the AGN feedback, which potentially has larger influence on 
the tSZ effect, although we do not expect it to significantly 
modify r{k,z). We shall investigate it in the future and quan- 
tify its possible impact on the tSZ tomography. 

There's also a possible limitation that may degrade tSZ 
tomography method using galaxies surveys. In our non- 
adiabatic simulation, extra gas dynamics may influence the 
galaxy-matter relation, resulting in deviations and stochastic - 
ities in the correlation coefficient r(k,z). Though we can't fig- 
ure out the exact influence of gas dynamics and galaxy selec- 
tion function, given little knowledge of recipes of galaxy for- 
mation at present, we test the possible effects with a generic 
halo occupation distribution (HOD) model in the appendix. 
The results of this generic HOD model favor tSZ tomography 
methods using galaxy survey in some aspects, while disfa- 
vor it in other aspects. Thus we can't pin down the influence 
of gas dynamics and selection function in galaxy surveys, 



and further investigations are needed. However, an optimistic 
news is that we DO demonstrate the robustness of the tSZ 
tomography, alternatively using tSZ survey with dark matter 
instead of galaxies, while we can measure the clustering of 
dark matter to high accuracy on scales of tSZ interest, using 
lensing tomography. We thus can expect a useful application 
of tSZ tomography in the future surveys of tSZ effect. 
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APPENDIX 

THE INFLUENCE OF EXTRA GASTROPHYSICS TO GALAXY-MATTER RELATION 

We would like to address a crucial point here in the appendix. The tSZ tomography method employed in this paper is based on 
an assumption that r mg {k) = 1 on relevant scales, i.e. galaxies are fully correlated with dark matter, and the correlation coefficient 
r mg is not sensitive to gastrophysics. This assumption assures that there's a deterministic map between the galaxy and matter 
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distribution, and thus the correlation coefficient of pressure and galaxy is equal to that of pressure and dark matter. Galaxy survey 
is actually an agency which introduces dark matter distribution in the tSZ tomography method. Therefore, the uncertainties in 
galaxy clustering of different galaxy samples may introduce stochasticities in estimating the underlying dark matter clustering 
features. In our simulations, when we emerge extra gastrophysics, dark matter would be influenced little, while galaxies, on the 
other hand, suffer from these gastrophysical processes. Exact scheme of gastrophysics is needed to model the galaxy formation 
for both simulations and theories, while it's beyond our knowledge so far. We thus have no idea how to incorporate detailed 
gas dynamics in both observations a nd simulation s given little knowledg e of the underlying gastrophysics. Nevertheless, earlier 
studies of galaxy-matter correlation (Seljak 2000; Guzik & Seljak 2001) showed that r is approximate to unity up to k ~ 1-2 



h Mpc , by cross correlating dark matter with galaxy samples selected by different selection criteria. In an alternative point of 
view, if we can measure matter power spectrum, we can circumvent the galaxy-matter cross corr elation while directly correlate 
dark matter and pressure power spectrum when doing tSZ tomography. [Seljak & Warrenl (120041) propo sed to use faint galaxies 
as dark matter tracer, as they found the linear bias of these galaxies is nearly constant. iPenl ( 120041) suggested to obtain the 
3D dark matter clustering using the cross correlation between galaxies with distance information and projected weak lensing 
dark matter maps. More recently. [Baldauf et al.l d2009j) proposed that by combining galaxy-clustering and galaxy-galaxy lensing 
measurements, they can reconstruct the dark matter power spectrum at a few percent accuracy up to k ~ 1 h Mpc -1 . 

However, we would like to investigate this correlation coefficient r mg in a more general way, say, by a HOD model in stead of 
looking into our two specific hydrodynamic simulations. Following the routine of halo model in lCoorav & S heth (20Q2|), we can 
calculate the power spectrum of dark matter as well as of galaxies. Given the halo model and the halo occupation distribution 
model (HOD) for galaxies, we can predict the matter power spectrum P m (k) and galaxy power spectrum P g (k), as well as the 
cross power spectrum P mg (k). The cross correlation coefficient is thus simply derived by r m Jk) = P mP (k)/ y/Pm(k)PJ k). We won't 
describe the details of how to calculate the power spectra, for which readers can refer to ICooray & Sheth! (120021) for a review. 
However, we present in the following paragraphs how we populate galaxies in halos with a generic HOD model. 

Though it's hard for one to exactly predict the galaxy distribution either by theories or by simulations, detailed semi-analytic 
galaxy formation models provides one with an experienced way to describe the galaxy distribution in the existence of generic 
gastrophysics. In this paper, we are not going to address various gastrophysics or thereby quantify the relation between galaxies 
and dark matter. As a simpl e approximation, we illu strat e the mean gal axy number (iVgai \M) in a halo given mass M by taking 
the parameters suggested by lSheth & Diaferiol (120011) and lCoorayl (12004) from the semi-analytic results, as follows 

IN \M\-I No for M min < M < M B 

{N s a \M) - | mM / MBr for M >M B ( A1) 

Here the parameters are M min = l0 n h~ 1 M G ,M B = 4 x 1O 12 /! _1 M , JVo = 0.7 and a = 0.8. We take this set of parameters as 
our fiducial model. Besides the mean occupation number of galaxies, we need to know as well the information of the second 
momentum of the galaxies. The u sually adopted Poisson probability distribution of galaxie s deviates from the resul ts from 
numerical data dBerlind et al.ll2003l) . we thus instead use the binomial distribution suggested by lScoccimarro et"aj] (1200 ll) 

(WgalCJVgal- D|M} 1/2 = /?(M) (JVgal|M) , (A2) 



here/3(M) = log s/M /M min M Q for M < Mq with Mq = 10 13 /z l M Q , and /3(M) = 1 otherwise. In ou r fiducial model, we assume the 
galaxies in halos trace the dark matter distribution, which follows the NFW (iNavarro et al.ll 1997b profile 

P {r)= < I vfl / 5 ' (A3) 
(r/r s )(l + r/r s y 

where p s is the density at a characteristic radius r„ which can be related to the virial radius r v j r with concentration c = r^/r s . A 
useful approximation for concentration is given bv lBullocket"aII(l200l 



c,„(M) = 9 



M 



-0.13 



(A4) 



There will be great differences between the HOD parameters of galaxy samples selected by various selection criteria, e.g. 
luminosity, color, history and morphology. However, we take into account the selection criteria and extra gastrophysics by 
considering two generic cases, which are framed by changing the number density of galaxies and changing the morphology of 
galaxy distribution. More specifically in our HOD model, we introduce these two effects by varying the six parameters of HOD 
adopted here, i.e. M m ; n , Mb, Mq, No, a and the concentration of galaxies c g . In practice, we calculate r mg by changing one of the 
parameters while fixing the others. In the existence of star formation and SN feedback, since it's hard for gas to cool and condense 
into stars, we expect less galaxies in the non-adiabatic simulation than the adiabatic one. We set the parameters such that the 
average galaxy number is less than that in the fiducial HOD model. We th us raise each time one of the mass thresholds, M m ; n , 
Mb and Mq as twice as that in the fiducial mod el, while reduce No by a half. lYang et aD (120051) found c g = 1 /3c,„ by investigating 
galaxy surveys, and lNagai & Kravtsovl (120051) found similar results in hydrodynamical simulations. Inspired by these works, we 
set the galaxy concentration c g as 1/3 or 2/3 as that of dark matter c m . For a, we reduce it to 0.7 so as to reduce the average 
galaxy number to test the effect. 

We show the result of the cross correlation coefficient r mg in Fig. [7] We find that reducing iVo doesn't change the r at all, and 
changing Mq only influence r by around 2 percent, except a 5% drop on scales around k ~ 5 h Mpc -1 . A more extended galaxy 
distribution in halos, i.e. c g = 1 /3c,„ would induce a 10% suppression of r mg at k ~ 7 h Mpc -1 . However, the suppression is 
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FIG. 7. — The sensitivity of the dark matter-galaxy cross correlation coefficient r mg to different HOD parameters. We change one of the parameters once at a 
time while keep others fixed during the calculation. 2M m ; n , 2Mb arid reducing a to 0.7 greatly influence r mg . However, it's possible that they cancel each other 
partly and reduce the stochasticity. Encouragingly, changing Nq, Mq, c g would introduce several percent influence on r on scales of interest of tSZ effect. 

at most 5% on scales of interest of tSZ, say, k < 3 h Mpc -1 . c g = 2/3c m , as expected, only reduce r mg by at most 3% on small 
scales. Actually an amplitude of galaxy number No won't introduce any difference of galaxy clustering. While galaxy number is 
insensitive to Mo or irrelevant to c g , the one halo term of galaxy power spectrum, i.e. on small scales is indeed influenced. 

The changes in M m ; n , Mg and a, on the other hand, lead to great influence in r mg . 2M m i n , i.e. increasing the mass threshold 
of host halo of central galaxy, reduces the average galaxy number by nearly a half, from 1.93 x 10~ 2 /r 3 Mpc~ 3 in the fiducial 
model to 9.93 x 10~ 3 /i 3 Mpc~ 3 , while 2Mb only changes that by 3%. Interestingly, they induce opposite effects in r mg , both with 
amplitudes up to 15%-20% on scales of interest. Though a = 0.7 reduces the galaxy number by less than 2%, it suppresses r mg by 
10% around k ~ 1 -2 h Mpc -1 . Thus from this simple model, it's possible that the stochasticity between dark matter and galaxies 
is large, and therefore the effectiveness of tSZ tomography would degrade. However, we hope that the change in r mg would not 
be so severe. Firstly, the parameters we adopt here are extreme, while in the non-adiabatic simulation, cooling and feedback 
influence mostly in the central regions of halos, and it would tend to influence small scales power spectrum below cluster scales. 
Secondly, they l ead to opposite effects to r mg , it's thus possible that the effects partly cancel out each other. In earlier works, 
iGuzik & Seli aid (120011) showed that galaxy-dark matter correlation coefficients of different galaxy samples, which are selected 
by luminosity, colour and star formation rate, converge to unity on scales above 1 h Mpc" . Also, they found in the faintest and 
most abundant samples, r is unity up to k ~ 3 h Mpc -1 , which is also indicated by iSeliak & Warre"n1 ( 120041) where they found 
small halos trace dark matter with constant bias. This is encouraging, as in practice, we can select faint galaxy sample to do tSZ 
tomography. After all, the stochasticity of galaxy-matter relation is an interesting issue, and need to be further investigated. 



